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ABSTRACT 

Emission lines from metals offer one of the most promising ways to detect the elusive 
warm- hot intergalactic medium (WHIM; 10 5 K < T < 10 7 K), which is thought to 
contain a substantial fraction of the baryons in the low-redshift Universe. We present 
predictions for the soft X-ray line emission from the WHIM using a subset of cos- 
mological simulations from the Overwhelmingly Large Simulations (OWLS) project. 
We use the OWLS models to test the dependence of the predicted emission on a 
range of physical prescriptions, such as cosmology, gas cooling and feedback from star 
formation and accreting black holes. Provided that metal-line cooling is taken into 
account, the models give surprisingly similar results, indicating that the predictions 
are robust. Soft X-ray lines trace the hotter part of the WHIM (T > 10 6 K). Wc find 
that the OVIII 18.97 A is the strongest emission line, with a predicted maximum 
surface brightness of ~ 10 2 photon s _1 cm~ 2 sr _1 , but a number of other lines are 
only slightly weaker. All lines show a strong correlation between the intensity of the 
observed flux and the density and metallicity of the gas responsible for the emission. 
On the other hand, the potentially detectable emission consistently corresponds to the 
temperature at which the emissivity of the electronic transition peaks. The emission 
traces neither the baryonic nor the metal mass. In particular, the emission that is 
potentially detectable with proposed missions, traces overdense (p > 10 2 p moa n) and 
metal-rich (Z > 10 _1 Zq) gas in and around galaxies and groups. While soft X-ray 
line emission is therefore not a promising route to close the baryon budget, it does offer 
the exciting possibility to image the gas accreting onto and flowing out of galaxies. 

Key words: method: numerical - intergalactic medium - diffuse radiation - radiation 
mechanisms: thermal - cosmology: theory - galaxies: formation 



1 INTRODUCTION 

Observations of the low-redshift Universe are unable 
to account for a large fraction of the baryons pre- 
dicted by cosmic microwave back ground (CMB) and nu- 
cleosynthesis measur e ments (e.g. Persic fe Saluccil 1 19921 : 
IFukugita et~ai] Il998l : IFukugita fc Peebles! l2004h . Cosmo- 
logical simulations predict that most of these elu- 
sive baryons reside in a diffuse, warm-hot intergalactic 
medium (WHIM, hereaft er) with temperatures i n the range 
10 s K < T < 10 7 K |Cen fc Ostrfkeil 1 19991 : iDave et all 
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2001 ; ICen fc Ostrikeil l2006t see iBertone. Schave fc Dolael 
20081 for a recent review). Theoretical models agree 



with numerical simulations that the relatively high tem- 
perature of the WHIM is produc e d by gravit a tional 



perature or the WHIM is produc e d by gra vitational 
heating dSunvaev fc Zel'dovicrJ Il972l: iNath fc Silkl l200ll : 



iFurlanetto fc Loebl |2004 iRasera fc Tevssieij l2006h . When 
large-scale structures form by gravitational collapse, cosmo- 
logical shock waves are driven into the intergalactic medium 
(IGM) and heat the gas to temperatures much greater than 
10 4 K. 



The detection of the WHIM in the low-redshift 
Universe requires observations in th e ultraviolet ( UV) 



and in the soft X-ray bands (see e.g. iBregmanl [20071 and 
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iProchaska fc TumlinsorJ |200S| for reviews). Unfortunately, 
UV and X-ray absorption lines are more technically chal- 
lenging to detect than optical lines, and require satellites 
outside the Earth's atmosphere. In addition, the low density 
of the WHIM implies that its surface brightness is extremely 
low, making it nearly impossible to detect in emission with 
current instruments such as Chandra and XMM-Newton 
and the Far Ultraviolet Spectroscopic Explorer FUSE 
JPierre et all |2000|; IKravtsov etldl |2002|; lYoshikawa et all 




et alJ 120051: lYoshida et al 120051 : lYosh ikawa & Sasaki 
20061; lUrsino fc Galeazzil 2006 ) . However. IWerner et al 



iFurlanetto et alJ 12004 ; lYoshikawa et alJ 12004 
Yoshida et all 120051 



claim to have detected continuum emission with 
XMM-Newton from a filament extending between the two 
large clusters Abell 222 and Abell 223 at 5<r significance, 
thanks to the orientation o f the filame nt along the line 
of sig ht. IZappacosta et al] (|2002l ) and IZappacosta et al] 
(2005) also found evidence for diffuse thermal emission 
from WHIM gas in two ROSAT field s, corresponding to an 
overdensity of galaxies at z ~ 0.4 (jMannucci et al.ll2007T ) 
and to the Scultor Wall. 

Absorption lines are less biased towards high gas 
densities and may therefore be more sensitive to the 
low-density WHIM. There are a few controversial claims 
of detection of O VII abso r ption lines in the spectrum 
of Mrk 421 llNicastro et all 120051 : but see iKaastra et al] 
l200d: iRasmussen et all 120071 ) and from the Sculptor Wall 
|Buote et al, ~ l2009l ). Future satellites such as the Inter- 
national X-ray Obser vator ^l (IXO hereaf ter, see also 
iMarkevitch et al.) 120091 and iBregmanl I2007T ) and Xeniffl 
jHartmann et al. 2009) will be ideal to detect a large num- 
ber of WHIM absorption lines in the spectra of QSO and 
gamma-ray bursts, respectively (e.g . Hcllstcn et al.l 1998; 



Perna fc Lo"eblll998l; IKravtsov et aT.ll20q2l;IChen et al.l | 2003 | 
Cen fc Fandl200d:lKawahara et al.ll2006l : IPaerels et al-lbool : 
Branchini et al] 12009 ) 



Although X-ray absorption lines are a powerful tool to 
investigate the physical properties of the WHIM, the infor- 
mation they yield has the limitation of being 1-dimensional, 
unless close quasar pairs can be found in large numbers. 
While absorption lines are not ideal to reconstruct the 3- 
dimensional distribution of the gas, emission lines do just 
that: they effectively map the 2-dimensional distribution of 
the WHIM on the sky, and by means of their spectral red- 
shift they trace the distribution of the gas along the 3rd 
dimension. Detecting emission from the WHIM is likely to 
be much more challenging than detecting absorption, but 
the reward is comparatively large. It is therefore paramount 
to predict accurately what level of emission can be expected, 
in order to build instruments and satellites adequate to the 
task. 

The spatial distribution of metals in the IGM, as traced 
by the WHIM emission, also provides important clues for un- 
derstanding the metal enrichment history of the Universe. 
The distance travelled by metals between the sites of star 
formation and the lower density IGM can constrain the en- 
ergy budget of processes such as galactic winds and AGN 
feedback. 



1 http://ixo.gsfc.nasa.gov/ 

2 http://sms.msfc.nasa.gov/xenia/ 



Table 1. List of the emission lines considered in this work. The 
first column shows the name of the emission line, the second and 
the third columns the wavelength and the energy of the transition, 
respectively. The last column identifies the type of ion responsi- 
ble for the transition (e.g. "H" means H-like) and the letter in 
parentheses indicates the transition for the case of triplets from 
He- like atoms: resonance (r), forbidden (f) and intercombination 



Ion 


A (A) 


E (eV) 


Descript 


CV 


40.27 


307.88 


He (r) 


CVI 


33.74 


367.47 


H 


N VI 


29.53 


419.86 


He (r) 


NVII 


24.78 


500.24 


H 


OVII 


21.60 


573.95 


He (r) 


OVII 


21.81 


568.74 


He (i) 


OVII 


22.10 


560.98 


He (f) 


OVIII 


18.97 


653.55 


H 


NelX 


13.45 


921.95 


He (r) 


NeX 


12.14 


1022.0 


H 


MgXII 


8.422 


1471.8 


H 


Si XIII 


6.648 


1864.9 


He (r) 


SXV 


5.039 


2460.5 


He (r) 


FeXVII 


17.05 


726.97 


Ne 



In this work, we use a subset of runs taken from 
the Over -Whelmingly Large Simulations (OWLS, hereafter) 
project (|Schave et al] 120101 ) to investigate which emission 
lines in the soft X-ray band best trace the low-redshift 
WHIM emission. We will do the same for UV emission in a 
companion paper (Bertone et al., in preparation). We esti- 
mate the level of the emission in a simulated cubic region 
of 100 h^ 1 comoving Mpc on a side and discuss the observ- 
ability of the WHIM by future telescopes such as IXO, the 
Explorer of D iffuse emission a nd Gamma-ray burst Explo- 
sions (EDGE. \p\ro et alj|2009h and Xen ia. 

The OWLS runs (|Schave et"ai1l2010h are well-suited for 
this study and represent a significant step forward with re- 
spect to previous studies. OWLS couples a large simulated 
dynamical range with a large set of models in which the 
implementation of many uncertain physical processes have 
been varied. This special feature allows us to investigate the 
dependence of the results on different physical processes. In 
this paper we will focus on the strongest soft X-ray emis- 
sion lines we have identified, which are listed in Table [I] 
the most important being OVIII A 18.97 A and OVII A 
21.60 A. When more than one line corresponds to a given 
transition, we consider only the strongest line. 

This paper is organised as follows. We briefly describe 
the OWLS reference model and the method to calculate the 
line emission of the WHIM gas in Section [2] Section [3] con- 
tains predictions for various emission lines for our reference 
model. Sections U and [5] are devoted to understanding the 
effects of varying the angular resolution of the detector and 
the redshift dependence of the flux intensity, respectively. 
Section HJ] investigates what kind of gas produces most of 
the emission and Section [7] studies the dependence of the 
results on the physical model. We discuss the prospects for 
observing soft X-ray emission lines with future space-borne 
telescopes in Section [8] and we summarize our conclusions in 
Section[9] Finally, we demonstrate in the Appendix that our 
conclusions are robust with respect to changes in both the 
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Table 2. Adopted solar abundanc es, f rom I AU ende P rieto et all 
j200ll l, I AUende Prieto et alj J2002T) and lHolweeej feOOll) . 



Element 


m /ng. 




Element 


rii/nu 


H 


1 




Mg 


3.47xl0~ 5 


He 


0.1 




Si 


3.47xl0~ 5 


C 


2.46x10" 


4 


S 


1.86X10- 5 


N 


8.51x10" 


-5 


Ca 


2.29xl0~ 6 





4.90x10" 


4 


Fe 


2.82X10"" 4 


Ne 


1.00x10" 


4 







numerical resolution and the size of the simulation volume. 
We also show in the Appendix how the results vary with the 
band width. 

Unless stated otherwise, we will assume the Wilkinson 
Microwa ve Anisotropy Prob e year 3 (WMAP3) ACDM cos- 
mology (jSpergel et al.l [20071 ) with parameters fi m = 0.238, 
n b = 0.0418, fi A = 0.762, n = 0.951, and <r 8 = 0.74. The 
Hubble constant is parameterised as Ho — 100 ft. - kms - 
Mpc - , with h — 0.73. The adopted solar abundance is 
Zq = 0.0127, corresponding to the value obtained using the 
default abun dance set of CLOUD Y (version c07.02.02, last 
described in iFerland et alj Il998h . Thi s abundance set (see 
Table [2t combines the abun d ances of lAllende Prieto et al.l 
|200lOA"Iiende Prieto et al.1 |2002l ) and Iflolwegerl (120011 ) 
and assumes UHe/fiH = 0.1. We note here that the abun- 
dances adopted by CLOU DY may differ b y significant factors 
from those estimated bv lLodderj J2003TI. In particular, the 
abundance of iron given bv lLoddera (|2003h is almost a factor 
of ten smaller than in CLOUDY, that of oxygen about 20 per 
cent higher. This should be kept in mind when comparing 
predictions from different studies. We stress, however, that 
the assumed solar abundances play no role when we use the 
abundances predicted by our simulations. 



2 THE SIMULATIONS 



The OWLS project (|Schave et alj|20ld) consists of a suite of 
more than fifty large, cosmological, hydrodynamical simula- 
tions that have been produced using a significantly extended 
version of th e parallel PMT ree-SPH code GADGET III (last 
described in ISpringeil 120051 1. In addition to hydrodynamic 
forces, the simulations include new prescriptions for radia- 
tive cooling, star formation, chemodynamics, and supernova 
feedback. Feedback from AGN is included in only one of the 
models (see Section [7]). Sections 0-[H] will only make use of 
the OWLS reference model, but we will investigate a range 
of other models in Section [7] 

The sub-set of the OWLS runs used in this work have 
boxes of size 100 comoving /t _1 Mpc, and contain 512 3 
particles of both gas and dark matter. Comoving gravita- 
tional softenings are set to 1/25 of the mean comoving inter- 
particle separation, but are limited to a maximum physical 
scale of 2 kpc (which is reached at z = 2.91). The 
initial gas particle mass is 8.66 x 10 7 h' 1 Mq. The mi- 
tial conditions were gener ated using cmbfast (version 4.1; 
ISeliak fe Zaldax riagalll996f l and evolved to z = 127 using the 
Zeldovich approximation from an initial glass-like state. 

Radiative cooling is implemented following the pre- 



scriptions of IWiersma et all ( 2009a Jf[ In brief, net radia- 
tive cooling rates are computed element-by-element in the 
pre sence of the cosmic micr owave background (CMB) and 
the lHaardt fe Madaul ([200 il l model for the UV/X-ray back- 
ground radiation from quasars and galaxies. The contribu- 
tions of the eleven elements hydrogen, helium, carbon, ni- 
trogen, oxygen, neon, magnesium, silicon, sulphur, calcium, 
and iron are interpolated as a function of density, tem- 
perature, and redshift from tables that have been precom- 
puted using the publicly ava ilable photo-ionizati on package 
CLOUDY (last described bv lFerland et al J [19981. assuming 
the gas to be optically thin and in (photo-)ionization equi- 
librium. l^i£_shi n J^ions model reionization by 'switching 
on' the lHaardt fe Madaul (|200ll l metagalactic UV/X-ray- 
background at z — 9. This has the effect of rapidly heating 
all the gas to temperatures of ~ 10 4 K. 

Star formation is tre a ted fo llowing the prescription 
of ISchave fe Dalla Vecchial j2008). Gas with densities ex- 
ceeding the critical density for the onset of the thermo- 
gravitational instability (nrr ~ 10 -2 — IP - 1 cm" 3 ) is ex - 
pected to be multiphase and star-forming (|Schavd |2004| 1. 
An effective equation of state (EOS) is imposed with pres- 
sure P oc p 7off for densities iih > 1h = 0.1 cm~ 3 , nor- 
malised to P/k = 1.08 x 10 3 cm- 3 K at the threshold. We 
use 7 c ff = 4/3 for which both the Jeans mass and the ra- 
tio of the Jeans length to the SPH kernel are independent 
of the density, thus preventing spu rious fragmentation due 
to a l ack of numerical resolution (ISchave fe Dalla Vecchial 
12008( 1 . ISchave fe Dalla Vecchial (|2008l l showed how the ob- 
served Kennicutt- Schmidt law, E* = A(E g /l Mope -2 )™, 
with E* the star formation rate surface density and E g the 
gas surface density, can be analytically converted into a pres- 
sure law, which can be implemented directly into the simu- 
lations. This has the advantage that the parameters are ob- 
servables and that the simulations reproduce the input star 
formation l aw irrespective of the assumed equation of state. 
We use the ISchave fe Dalla Vecchial (|2008l l meth od, setting 
A = 1.515 x 10~ 4 M©yr _1 kpc~ 2 and n = 1.4 (|Kennicuttl 
1998, note that we renormalized the observed relation for a 
Chabrier IMF). 

As described in W iersma et al.l (|2009bl l , the simulations 
follow the timed release of all 11 different elements that con- 
tribute significantly to the cooling rates by massive stars 
(Type II supernovae and stellar winds) and intermediate 
mass stars (Type la supernovae an d asymptotic gi ant branch 
stars), assuming a Chabrier IMF |Chabrierll2003l l spanning 
the range 0.1 to 100 Mq. 

Energy injection due to supernovae is included 

through kinetic feedba ck following the prescription of 

iDalla Vecchia fe Schavel |2008l l. Core-collapse supernovae 
locally inject kinetic energy and kick gas particles into 
winds. The feedback is specified by two parameters: the ini- 
tial mass- loading 77 = m w /m*, which describes the initial 
amount of gas put into the wind, TTI-uj , JUS a function of the 
local SFR m„, and the wind velocity d w . We use r\ = 2 and 
w w = 600 kms" 1 , which corresponds to 40% of the total 
amount of supernova energy being available in the form of 



3 Using their equation (3) rather than (4) and CLOUDY version 
05.07 rather than 07.02. 
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kinetic energy and yields a cosmic star formation history i n 
good agreement with the observations (ISchave et al.ll2010h . 

Finally, we note that W lersma et all (l2009bF showed 
that the mass-weighted mean metallicity of the WHIM is 
about 0.1Z@ for the reference model. 



2.1 Emissivity tables 

The intensity of the emission lines is calculated by inter- 
polating tables of the line emissivity as a function of gas 
temperature, hydrogen number density and redshift. The 
calculation of the emission of individual particles in the sim- 
ulation is described in Section 12.21 The tables have been 
created using CLOU DY (version c07.02.02, last described in 
iFerland et al 1 ll99Sl and released in July 2008) and include 
about 2000 emission lines for the 11 elements tracked explic- 
itly by our simulations. CLOUDY was run using the same 
assumptions as were used for the calculation of the cooling 
rates: optically thin gas in (photo- )ionisation equilibr i um in 
the presence of the CMB and the lHaardt fc Madaul (fcOOlT l 
metagalactic UV/X-ray background. The CLOUDY tables 
were computed assuming solar abundances (see Table [2]), 
but are scaled to the abundances predicted by the simula- 
tion (see Section T2.2[l , The tables have been created for the 
temperature range 10 2 K < T < 10 8 ' 5 K and the hydrogen 



number density range 10" s cm -3 < n H < 10 cm 



The 



temperature is sampled in bins of ALog 10 T = 0.05 and the 
hydrogen number density in bins of ALog 10 7iH = 0.2. 

The emissivities of a selection of the strongest individ- 
ual lines, listed in Table [1] are shown in Fig. [T] as a func- 
tion of temperature, for constant hydrogen number densi- 
ties riH = 10 -6 , 10 -3 and 1 cm -3 and assuming solar abun- 
dances. More lines are included in our tables, but since their 
intensities are much weaker than for those listed, we do not 
consider them for this work. However, these weaker lines are 
crucial for creating synthetic spectra of the WHIM emission, 
because their cumulative contributions can significantly in- 
crease the level of the continuum. 

Collisional ionisation is the dominant ionisation pro- 
cess at the highest densities and temperatures, while the 
contribution of photo-ionisation is largest at low den- 
sity and at low temperatures. The wings of the emis- 
sivity curves at T < 10 6 K and tih < 10~ 3 cm -3 are 
signatures of photo-ionisation. The assumption of ionisa- 
tion equilibrium, also used to calculate the cooling ta- 
bles |w lersma et all l2009al l . is well justified for photo- 
ioni sed regions and dense gas in the centres of clusters 
(see iBertone. Schave fc D olag 2008 for a review). However, 
non-equilibrium processes may become important in the 
WHIM and in the outer regions of groups and clusters 
dYoshida et al.ll2005l; lYoshikawa fc SasakilliobHilCen fc Fand 
l2006l ; lGnat fc Sternberdl2007h . lGnat fc Sternberd (|2009h .~in 
particular, have recently shown that non-equilibrium pro- 
cesses affect the cooling of shock-heated gas with temper- 
atures with 10 4 K < T < 10 7 K, and that the amplitude 
of the deviations from ionisation equili brium increases with 
the ga s metallicity. On the other hand, lYoshikawa fc Sasakil 
(2006) found the effect of non-equilibrium ionisation to be 
small for the subject of this paper: metal-line emission from 
the WHIM. 

The emissivitie s show n in Fig. [1] are comparab le to those 
of lYoshikawa et all ((2003) and lFang et al. I (|2005l ). with the 



exception of th e FeXVII l i ne, wh ich is significantly weaker 
than found by iFang et al.l (120051 ). In general, we find that 
CLOUDY predicts a Fe XVII emissivity a factor of about two 
small er than t he Raymond-Smi th code l|Ravmond fc Smith! 
ll9T7f) used bv lFang etaH (120051 ). Part of the reason for this 
discrepancy could be that two of the five strongest Fe XVII 
emission lines are missing from the CLOUDY line list, namely 
A15.01 A and A17.10 A. This might be an issue for the cool- 
ing of metal-rich gas at T ~ 10 7 K, where the FeXVII line 
can contribute up to 30 per cent of the total cooling rate. 
On the other hand, this lack of iron lines might be balanced 
by the high iron abu ndance as s umed in CLOUDY, compared 
to that estimated bv lLodders! (|2003l ). 



2.2 Line emission calculations 

We investigate the WHIM emission in the soft X-ray band 
by creating synthetic maps of the surface brightness of in- 
dividual emission lines. We only include emission from gas 
with densities tih < 0.1cm -3 , because our simulations lack 
the resolution and the physics to model the emission from 
higher density g as, which is expected to be mu ltiphase (see 
the discussion in iDalla Vecchia fc Schavell2008l ). 

We calculate the luminosity of emission line I for each 
gas particle i as a function of its gas mass m gas ,i, density p;, 
hydrogen number density nu,i and element abundance X y ^ 



L l} i (z) = Si, !j0 (2,Ti,n H ) 



Pi X y Q 



(erg s x ) , (1) 



where e^,© (z, Ti, tih) is the emissivity of the line I for solar 
abundances, in units of erg cm" 3 s -1 , bi-linearly interpo- 
lated (in logarithmic space) from the CLOUDY tables de- 
scribed in Section 12.11 as a function of the particle tem- 
perature Log 10 Ti and hydrogen number density Log 10 riH,i 
at the desired redshift. X Vl i is the mass fraction of ele- 
ment y and X y Q its solar value. We use the "smoothed " 
element abundances described in W lersma et al.1 (|2009bl ). 
which are defined as the ratio of the SPH smoothed mass 
density of the element and the SPH smoothed total mass 
density, Z y> i = p y ,i/ 'pi. While the use of smoothed abun- 
dances, which were also used to calculate the cooling rates 
during the simulations, reduces the effects of the lack of 
metal m ixing inherent to SPH , it does not solve the prob- 
lem (see lWiersma et al.ll2009bl for a discussion). 

The corresponding particle flux Fij is given by 



Fi,i = 



(1 + z) (photon s 1 cm 



(2) 



with Dl the luminosity distance, hp the Planck constant, 
c the speed of light and A; the rest-frame wavelength of 
the emission line I. Finally, we project the fluxes on to a 
2-dimensional grid whose number of pixels is determined 
by the angular resolution we require and by the redshift we 
want to "observe" . We will investigate the effect of changing 
the pixel size in Section [4] The projection is done using a 
flux-conserving SPH interpolation scheme. 

Surface brightness (Sb, in all figures) maps are com- 
puted by dividing the flux in each pixel by its solid angle 
f^pixci = 2-7T [1 — cos (#rad/2)], where i? ra d is the pixel size in 
units of radians. At z — 0.25 our simulation box corresponds 
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Figure 1 . The normalised emissivity, e / n| , in units of erg cm 3 s — 1 , of a selection of soft X-ray emission lines as a function of temperature, 
for constant hydrogen number density and assuming solar abundances in the presence of the redshift z = 0.27 UV/X-ray background. 
The left, middle, and right panels show results for njj = 1, 10 — 3 and 10 — 6 cm - 3 , respectively. If the gas is collisionally ionised, the 
normalised emissivity is independent of the gas density. Cool diffuse gas is mostly photoionised, while collisional ionisation equilibrium 
dominates for the highest temperatures and densities. 



to a field of view of about 8 degrees and an angular resolu- 
tion of 15" requires a grid with 1926 2 pixels. At z — 0.25 this 
corresponds to a physical scale of about 42 h~ kpc, which is 
about 21 times larger than the gravitational softening scale 
used in the simulation. Unless stated otherwise, the emis- 
sion maps shown were created with a pixel size correspond- 
ing to an angular resolution of 15" on the sky. We divide 
each simulated box of 100 Mpc into 5 slices that are 
each 20 Mpc thick and create as many emission maps, 
which we combine to calculate the flux probability distribu- 
tion function (PDF, hereafter) corresponding to this angular 
resolution and slice thickness. We neglect the effect of pecu- 
liar velocities when creating the slices, and select particles 
only according to their position within the simulation box. 
At z = 0.25, a slice thickness of 20 h' 1 Mpc corresponds to 
an energy resolution 

^ = ^ = ^M-^-«6xlO-, (3) 

E c c 1 + z w 

which is comparable to the spectral resolution proposed for 
IXO and EDGE. In Appendix I A3I we demonstrate that the 
flux probability density function is proportional to the slice 
thickness, at least for AL > 5 h -1 Mpc and sufficiently large 
fluxes (greater than 10 -5 photon s _1 cm -2 sr _1 for OVIII 



and our fiducial angular resolution of 15"). This is because 
the probability of a pixel having a certain flux is propor- 
tional to the volume probed, provided the slice is sufficiently 
thick to encompass the emitting objects and sufficiently thin 
for projection effects to be unimportant (this last condition 
breaks down at very low flux levels). 



3 EMISSION MAPS 

In this Section we present emission maps for the strongest 
soft X-ray emission lines we have identified as WHIM trac- 
ers. All PDFs are calculated using the full simulation box 
with size 100 h^ 1 Mpc. However, for clarity we show emis- 
sion maps for only a fraction of the full box. In particular, 
we focus on a smaller area of 14 h' 1 Mpc on a side that 
contains one of the largest groups in the simulation, with a 
total mass of about 7 x 10 14 Mq, a smaller group and a few 
dense filaments. This choice is motivated by the fact that 
the strongest emission comes from high-density regions such 
as groups and filaments. 

To illustrate the morphology of the emission, Fig. [2] 
shows the O VIII emission from a slice of the full simulation 
box in the left panel and two consecutive zooms into the core 
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Figure 2. Zoom into an O VIII 18.97 A emission map of a high-density region at z = 0.25. The object at the centre of the displayed 
region is one of the largest haloes in the simulation, with a total mass of a few times 10 14 Mq. The slice thickness is 20 h _1 comoving 
Mpc (which corresponds to 1772 kms -1 ). From left- to-right, the pixel sizes are 100", 10" and 1" and the angular sizes of the regions 
shown are 8 degrees, 48' and 4.8' (corresponding to 100, 10, and 1 /i" 1 comoving Mpc), respectively. The blobs in the right panel are 
resolved substructures. 



of a massive group in the companion panels. The displayed 
regions decrease in size by factors of ten, ranging from 100 
/t" 1 Mpc to 1 h~ x Mpc. The strongest O VIII emission is 
concentrated in the high-density knots at the intersections 
of filaments, where groups reside, while the emission from 
cooler and more diffuse gas in filaments and voids is several 
orders of magnitude weaker. For comparison, the intensity 
of the unre solved X-ray background is ~ 10~ 2 photon s" 1 
cm" 2 sr -1 l|Hickox fc Markev itch 2006, see also Section [8]) 
and roughly corresponds to orange in the colour scale. As we 
will see in the following, the surface brightness of the WHIM 
is higher than the background only in relatively overdense 
regions. 

Fig. [3] shows emission maps for a selection of twelve of 
the strongest emission lines in the soft X-ray band, listed in 
Table [T] The surface brightness PDFs for the same lines are 
shown in Fig. [4] in units of photon s _1 cm -2 sr -1 , and in 
Fig. [5] in units of erg s" 1 cm" 2 sr" 1 . To convert a flux F p h 
in units of photon s _1 cm" 2 sr" 1 to a flux F erg in units of 
erg cm" 2 s" 1 , we use the conversion 



F C1 



F 



P h 



A (1 + z) ' 



(4) 



with Ao the rest-frame wavelength of the transition (e.g. 
18.97 A for O VIII emission). 

Here we consider only the strongest lines from 
Hydrogen-like atoms and the resonance lines in triplets from 
Helium-like atoms, and we defer the discussion of the be- 
haviour of lines within the triplet to Section T3. II These lines 
have been chosen on the basis of their strength and de- 
tectability in the soft X-ray band, but are by no means the 
only important lines. Although we only show results for one 
iron line from the Fe XVII ion, we note that the soft X-ray 
band is a very rich band for iron lines from a multiplicity of 
other higher-order ions, including, e.g., Fe XVIII and Fe XX. 



Most lines from these ions trace slightly hotter gas than the 
FeXVII lines, but nevertheless represent an extremely im- 
portant cooling channel for hot gas in dense environments. 
Lines from Fe XVI are instead emitted at longer wavelengths 
and fall outside the soft X-ray band. 

Figs. 0(4] and [5] clearly indicate that O VIII A18.97 A is 
the strongest emission line and that for all lines the strongest 
fluxes are concentrated in very small regions or, equivalently, 
in a small fraction of the pixels in the maps. All lines are 
weaker than OVIII by a factor of a few or more. This is a 
consequence of the OVIII A18.97 A line emissivity, which 
is highest of all lines, and of the high relative abundance 
of oxygen. The emissivity of O VIII peaks at a temperature 
T ~ 3 x 10 6 K, and remains higher than for any other line at 
higher temperatures. The line emissivities shown in Fig. [1] 
give a good indication of what lines dominate the emission. 
However, the total emission in each line is not only deter- 
mined by the emissivity itself, but also by the fraction of the 
gas mass and metallicity in the temperature-density plane 
that corresponds to each emissivity value. 

Fig. [3] illustrates how some lines, namely C V, C VI, 
O VII and Ne IX, preferentially trace the cooler gas in 
filaments and in the outskirts of groups, while others, 
namely NeX, MgXII, Si XIII and SXV better trace the 
hot medium at the centres of groups. This is a direct con- 
sequence of the fact that the peak of the emissivity shifts 
to higher temperatures with increasing atomic number and 
ionisation state, as can be seen in Fig.[T] The C V, O VII and 
Ne IX lines, whose emissivities peak at about 10 6 K, 2x 10 6 K 
and 4 x 10 6 K, respectively, are best-suited to investigate the 
diffuse gas in filaments and in the outskirts of groups. How- 
ever, their emission declines drastically in the centres of large 
haloes, where the gas is very hot and more highly ionised 
atoms dominate (see e.g. the centres of the groups shown in 
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Figure 3. Maps of the emission for a sample of twelve different soft X-ray lines, listed in table [T] at z = 0.25. The figure shows a zoom 
of a region, with a comoving size of 14 h _1 Mpc (1.12 degrees on the sky), that includes the largest group in the simulation. The density, 
temperature, and metallicity in this region are shown in Fig. 1151 All panels assume the same angular resolution (15", which corresponds 
to a physical size of 42 h~ 1 kpc) and slice thickness (20 h^ 1 comoving Mpc, which corresponds to 1772 kms -1 ). The maximum of the 
colour scale has been set to a fraction of the real maximum to enhance the emission of the weaker lines and of low-density regions. The 
O VIII line is the strongest emission line in the sample. Carbon, nitrogen, oxygen and neon lines best trace the distribution of the WHIM 
in filaments and in the outskirts of groups, while lines from elements with higher atomic numbers such as MgXII, Si XIII, SXV and 
FeXVII mostly trace hot and dense gas in massive haloes. 



Fig. [3} . Emission lines such as CVI, O VIII, NeX, MgXII, 
Si XIII and S XV are the best probes of the dense gas in 
groups and clusters, but do not trace the low density gas in 
filaments and the field as well as C V, O VII and NelX do, 
because their emissivities peak at higher temperatures. 

The Fe XVII line is an exception to this trend because, 
even though its emissivity peaks at T ~ 6 x 10 6 K, its value 
is only high in a narrow temperature range. In addition, the 
spatial abundance of iron differs significantly from that of 
other elements because of the delay in the SN la enrichment. 
While iron is abundant in the intragroup medium and in the 



vicinity of galaxies, it is significantly underabundant in the 
diffuse IGM. The flux PDFs in Fig. Q] show that most lines 
produce fluxes higher than about 1 photon s _1 cm -2 sr -1 
in at least a few pixels. 

Fig. [6] shows the relative intensity of lines with respect 
to each other. When the contours lie above the dashed di- 
agonal line, the emission line on the t/-axis produces the 
highest flux, when below, the line on the a>axis is strongest. 
In general, the O VIII and C VI lines produce the highest 
fluxes, both among the lines that trace the hottest gas and 
among those that trace cooler gas. Lines from helium-like 
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Figure 4. The flux PDFs of the sample of soft X-ray emission lines that is also shown in Fig.[3]and that is listed in table[T] at z = 0.25. 
The left panel shows results for CV, CVI, NVI, NVII, OVII and OVIII, and the right panel for NelX, NeX, MgXII, Si XIII, SXV 
and FeXVII. The main plotting area shows only the high flux tail of the distributions, while the full PDFs are shown in the inset. The 
pixel size is 15" (which corresponds to a physical size of 42 h- 1 kpc) and PDFs are calculated using five slices through the simulation box 
that are each 20 h -1 comoving Mpc thick (which corresponds to 1772 kms" 1 ). In Appendix lA3l wc show that for sufficiently large fluxes 
(> 10 — 5 photon s -1 cm -2 sr _1 for OVIII) the PDF is proportional to the slice thickness. The OVIII line is the strongest emission 
line in the sample, followed by C VI and Ne X. The Fe XVII can reach high fluxes, but only in a limited number of pixels compared 
with other lines (see the inset of the right panel). This is due to the different spatial distribution of iron in the IGM and to the sharply 
declining values of the Fe XVII emissivity curve outside a narrow range around its peak temperature. The intensity of the unresolved 
X-ray background (~ 10 -2 photon s" 1 cm -2 sr — x , discussed in Section[8]l corresponds to the minimum value shown on the x-axis. 




Figure 5. As Fig. [4] but with the emission line flux expressed in erg s _1 cm -2 sr _1 , instead of photon s — 1 cm - 2 sr _1 . The relative 
strength of the lines varies slightly when converted to these units and is more easily comparable to the flux limits proposed for future 
missions such as IX O and EDGE. 
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Figure 6. Comparison of the intensities of the lines listed in Table [T] The contours show the number density of pixels and are logarith- 
mically spaced by 1.7 dex. The distributions have been calculated using five 20 hT x comoving Mpc slices through the simulation box 
at 2 = 0.25, assuming an angular resolution of 15". Both axes show the line surface brightness Log^Se in units of photon s — 1 cm -2 
sr -1 . The dashed diagonal line in each panel indicates equal fluxes. The O VIII and C VI lines produce the highest fluxes. Lines whose 
emissivities peak at similar temperatures (i.e. panels near the diagonal running from the bottom-left to the top-right of the plot) are 
strongly correlated. 



atoms produce higher fluxes than O VIII and C VI in a rel- 
atively small number of pixels. These correspond to regions 
of the IGM with relatively low temperatures, where helium- 
like atoms are the dominant species. The lines in Fig. [B] are 
sorted by atomic mass and then by ionisation state, which is 
nearly the same as sorting by the temperature for which the 
emissivity peaks (see Fig. [T]). This explains why the contours 
are much narrower for panels near the diagonal running from 



the bottom-left to the top-right of the figure than they are 
for panels near the top-left and bottom-right of the plot: 
panels near the diagonal correlate lines for which the emis- 
sivity peaks at similar temperatures. Fig. [5] illustrates the 
distribution of line ratios that can be expected in surveys 
of WHIM emission and may help to confirm the validity of 
line identifications. 

Although the actual difference in flux between emission 
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Figure 9. As Fig. g] but for the flux PDFs of the OVII triplet 
emission lines. The fluxes from the resonance and forbidden lines 
are about equal for Log 10 F > —3 photon s~ x cm -2 sr — x , but 
at lower fluxes the forbidden line becomes the dominant emission 
channel, while the intercombination line becomes as strong as the 
resonance line. 



lines can be as high as an order of magnitude, all emission 
lines in our set might be detectable in the centres of galaxy 
groups and in the outskirts of clusters. 



3.1 The O VII triplet 

In this Section we discuss the properties of line triplets 
from Helium-like atoms, focussing in particular on the O VII 
triplet (see Table [TJ, which should be the easiest to detect. 
The properties of the lines in triplets are important because 
the column density of a gas can be estimated using the ratio 
of the intensity of the resonance line of Helium-like atoms 
and other lines in the triplet. This might be possible, for ex- 
ample, when at least two out of the three lines are detected. 

The lines in the triplet correspond to the atomic tran- 
sition \s2p — > Is 2 and in particulaJJ 



(i) resonance line: x Pi — ► 1 So 

(ii) intercombination line: 3 P2,3 

(iii) forbidden line: 3 Si So 



So 



Of the three main emission lines of Helium-like atoms, 
the resonance line is usually of comparable strength to the 
forbidden line at high density, while the intercombination 
line is usually the weakest. The forbidden line, however, be- 
comes the dominant emission line in low density environ- 
ments, where gas is mostly photo-ionised and radiative in- 
tercombination is the dominant line emission process. 



4 Ralchenko Y., Kramida A.E., Reader J., and NIST ASD Team 
(2008). NIST Atomic Spectra Database (version 3.1.5). Available: 
http:/ /physics. nist.gov/asd3 National Institute of Standards and 
Technology, Gaithersburg, MD. 



Photons at the wavelength of the resonance line can 
be resonantly scattered in random directions when pass- 
ing through an intervening medium, while photons at the 
intercombination and forbidden frequencies have negligible 
cross-sections. As such, only the intensity of the resonance 
line can be attenuated along the line of sight, giving rise 
to absorpti on in the spectr a of ba ckground objects. Inter- 
estingly, as IChurazov et al.l (|200ll ) point out, the photons 
of the cosmic X-ray background are scattered in the same 
way, therefore boosting the line emission of the WHIM 
gas in filaments where the scattering occurs, compared to 
the thermal emission from the gas itself. This effect might 
improve the chances o f detecting the WHIM in emission 
l|Yoshikawa et alj|200j ). 

In Fig. [7] we show the flux predicted for the resonance, 
intercombination and forbidden lines in the O VII triplet. 
The flux PDF of the same lines is shown in Fig. [51 while 
Fig. [8] compares the flux intensity of pairs of lines on a pixel 
by pixel basis. The intensity of the intercombination line is 
about 2-10 times smaller than those of the forbidden and 
resonant lines. At the highest fluxes, F ^ 10~ 3 photon s _1 
cm~ 2 sr -1 , corresponding to high-density regions, as we will 
show in the next Sections, the intensity of the resonance line 
matches that of the forbidden line well, and the shapes of 
the PDFs are basically the same. At lower densities, the 
fluxes predicted for the forbidden line become higher than 
those of the resonance line by up to a factor of 10. This can 
be seen, for example, as a blue glow in voids in the right 
panel of Fig. [7J while the flux of the resonance line in the 
left panel is lower than the minimum of the colour scale. As 
a consequence, the ratio of the intensity of the resonance 
and the intercombination line, as well as the ratio of the 
resonance and the forbidden line, varies (left and middle 
panels of Fig. [8]), while only the ratio of the intercombination 
and forbidden lines is constant (right panel of Fig. [S| • 

In the absence of absorption by an intervening medium, 
both the forbidden and the resonance lines might be used 
to investigate the properties of WHIM gas. However, in low- 
density regions such as filaments and in the presence of an 
optically thick layer of gas, the forbidden line might be a bet- 
ter tracer of the WHIM. On the other hand, the resonance 
line is the only useful line for comparisons with absorption 
studies, such as those planned for IXO. 



4 ANGULAR RESOLUTION 

In this Section we investigate the effect of varying the an- 
gular resolution of the detector. 

Fig. [TO] compares the PDFs of the measured flux of 
OVIII emission assuming angular resolutions of $ = 5", 
15", 1', 4' and 10'. The distributions roughly converge for 
low fluxes, but there is no convergence at the highest fluxes 
and reducing the pixel size continues to increase the highest 
observed flux, at least down to the smallest scale considered 
here. This is because larger pixel sizes prevent the detec- 
tion of very strong emission coming from structures with an 
angular size smaller than that used to observe them. When 
observing diffuse gas on scales larger then galactic scales 
(such as the WHIM), an angular resolution of about 10-15" 
is sufficient at z ~ 0.25 and does not imply a significant loss 
of information. However, higher angular resolution might be 
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Figure 7. As Fig. [3] but comparing the resonance (left panel), intercombination (middle panel) and forbidden (right panel) emission lines 
of the O VII triplet. The resonance and forbidden lines have similar strengths in dense regions with strong emission, but the forbidden 
line becomes the main emission channel in lower density regions with fluxes lower than about 10 — 3 photon s — 1 cm -2 sr - 1 . 
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Figure 8. As Fig. [6] but for the line intensities in the OVII triplet: intercombination line vs. resonance (left panel), forbidden vs. 
resonance (middle panel) and forbidden vs. intercombination (right panel). 



desirable when observing gas with more clumpy spatial dis- 
tribution on smaller scales, such as outflowing gas in galactic 
winds. For reference, in the OWLS cosmology a pixel size of 
15" subtends a comoving (proper) size of about 52 (42) h~ 
kpc at 2 = 0.25, 98 (66) h' 1 kpc at z = 0.5 and 175 (87) 
kpc at z — 1. 



5 EVOLUTION 

In this Section we briefly discuss how the intensity of the 
emitted flux varies with cosmic time. 

As mentioned in the previous Section, as redshift in- 
creases, the physical size that corresponds to a pixel of a 
constant angular size increases steadily, at least up to z = 1. 
Moreover, the intrinsic emission depends on the tempera- 
ture, density and metallicity of the gas itself, as well as on 



the UV/X-ray background, all of which evolve. In Fig. llll we 
show the flux PDF of O VIII emission at redshifts z = 0.25, 
0.5 and 1. The overall shape of the flux PDF changes lit- 
tle with redshift, but the change is most noticeable at the 
highest fluxes, where the maximum flux decreases with in- 
creasing redshift. This is at least partly due to the fact that 
at constant angular resolution and over the redshift range 
considered here, pixels correspond to larger physical sizes at 
increasing redshifts and, as shown in Section [3] the maxi- 
mum flux decreases with increasing pixel size. 



6 WHAT TYPE OF GAS DOMINATES THE 
EMISSION? 

In this Section we investigate what type of gas contributes 
most of the soft X-ray emission. We begin by showing 
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Figure 10. As Fig. [4] but for OVIII emission as a function of 
the angular resolution of the emission maps. The pixel sizes are 
i9 = 5" , 15", 1', 4', and 10', corresponding to physical sizes of 
14 h- 1 kpc, 42 h- 1 kpc, 166 h~ x kpc, 665 kpc and 1.661 
h Mpc, respectively. Results broadly converge at high angular 
resolution (i? < 1'), except for the highest fluxes at the centres of 
haloes, which arise from gas in clouds with physical dimensions 
smaller than the pixel sizes considered here. There is no clear 
convergence at low angular resolutions > 1'). 

the emission-weighted gas distribution in the temperature- 
density plane for several lines in Section f6. II In Section [6. 2 1 
we then correlate the O VIII flux from individual particles 
with their density, temperature and metallicity. Finally, we 
calculate the O VIII emission from gas in different density 
and temperature ranges in Section 15751 

6.1 Temperature-density plane 

Fig. [T5] compares the mass-weighted and the metal mass- 
weighted temperature-density distributions of gas at z = 
0.25. The vertical line at nu — 0.1 cm~ 3 indicates the den- 
sity threshold above which star-forming particles are put 
on a power-law equation of state (see Section [2} and the 
horizontal lines at T = 10 5 K and T = 10 7 K delimit the 
temperature range that defines the WHIM. 

The mass-weighted and the metal mass-weighted dis- 
tributions in Fig. [12] clearly indicate that the bulk of the 
metals do not trace the bulk of the gas mass. A large frac- 
tion of the gas mass resides in the low-density IGM with 
T < 10° K, while metals are more common in the moder- 
ately dense WHIM (see also IWiersma et al. I l2009bl ). Also, 
a significant fraction of metals are found in the relatively 
dense gas which populates the branch of the distribution 
with riH > 10 -2 cm -3 and T < 10 5 K. This gas is accreting 
onto galaxies and may eventually end up in the interstellar 
medium (ISM). 

Fig. [13] presents the emission-weighted temperature- 
density distributions of the gas, where the weight is given 



Figure 11. As Fig. [4] but for OVIII emission as a function of 
redshift. Results are shown at z = 0.25, 0.5 and 1. Maps are cre- 
ated assuming constant angular resolution i9 = 15" , which implies 
the physical size of a pixel increases with redshift from 42 h~ 1 kpc 
at z = 0.25 to 88 h^ 1 kpc at z = 1. The flux PDFs are nearly the 
same, but diverge at the highest fluxes. In particular, the highest 
fluxes are found at the lowest redshift. This is probably mostly 
due to the varying physical size of the pixels. 
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Figure 12. The mass-weighted (left panel) and metal mass- 
weighted (right panel) temperature-density distribution of gas at 
z = 0.25. The large difference between the two distributions im- 
plies that the bulk of the metals do not trace the bulk of the 
gas mass. The contours are logarithmically spaced by 1.5 dex. 
The horizontal lines at constant temperature highlight the WHIM 
range, while the vertical line at nrr = 0.1 cm -3 indicates the limit 
above which we impose an equation of state onto star-forming par- 
ticles. These dense particles are excluded from the calculation of 
the emission. For reference, the cosmic mean density corresponds 
to n H ~ 3.7 X 10~ 7 cm -3 . 
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Figure 13. The emission-weighted temperature-density distributions of gas at z = 0.25. The distributions are weighted by the emission in 
different line for different panels, as indicated in the bottom-left corner of each panel. From left to right and from top to bottom, the lines 
are: C V, C VI, N VI, N VII, O VII, O VIII, NelX, NeX, MgXII, Si XIII, SXV and FeXVII. The contours are logarithmically spaced by 
1.5 dex and represent equal emission levels in all panels. The horizontal dashed lines at constant temperature highlight the WHIM range, 
while the vertical line at njj = 0.1 cm -3 indicates the limit above which we impose an effective equation of state onto star-forming gas. 
We do not include emission from this high-density gas. For reference, the cosmic mean density corresponds to njj ~ 3.7 X 10~ 7 cm~ 3 . 
Lines from atoms with low atomic numbers are better tracers of WHIM emisison, while lines from heavier elements and from hydrogen- 
like ions mostly trace hotter gas with temperatures in excess of 10 7 K. Fe XVII emission is only present in a very narrow temperature 
range, while C VI and O VIII lines trace gas in a wide temperature range. 



by the emission in a different line in each panel. The lines 
we consider here are C V, C VI, N VI, N VII, O VII, O VIII, 
NelX, NeX, MgXII, Si XIII, SXV and FeXVII, as listed 
in Table [1] Comparison of the mass- weighted, metal mass- 
weighted, and the emission-weighted temperature-density 
distributions demonstrate that neither the bulk of the mass 
nor the bulk of the metals are good tracers of the bulk 
of the soft X-ray emission from metal lines. The peak of 
the emission-weighted distributions of all lines is shifted to- 
wards higher densities, compared to the mass- and metal 
mass-weighted distributions. In fact, the emission of virtu- 
ally all highly ionised ions emitting in the X-ray band con- 
centrates in a region of the temperature-density plane with 
7iH > 10 _ cm - 3 and T > 10 B K. While the emission from 
hydrogen-like atoms is strong also for temperatures in ex- 
cess of 10 7 K, emission from helium-like atoms such as C V, 
NVI and O VII peaks at T < 10 7 K and declines steeply at 
higher temperatures. 

Fig. 1131 demonstrate that, among the ones shown here, 
the best lines to trace the WHIM in the temperature range 
10 6 K < T < 10 7 K are CVI, NVI, NVII, OVII, OVIII 
and Ne IX. All other ions tend to trace hotter and somewhat 
denser gas with T > 10 7 K and n H > 10" 4 cm -3 . 



6.2 Correlations between luminosity and gas 
properties 

Fig. [T3] shows what densities and temperatures dominate 
the global emissivity in various emission lines, but much of 
this emission may come out at fluxes that are too faint to 
be detectable. It is therefore also of interest to investigate 



physical conditions in the gas that produces the brightest 
emission. 

Fig. [14] investigates quantitatively the dependence of 
CVI, OVII, OVIII and NeX emission on the gas prop- 
erties for a sample of emission lines by showing the median 
values of the particle density (left panel), temperature (mid- 
dle panel) and metallicity (right panel) that produce a given 
particle luminosity. Since we consider individual particles for 
this plot, the density, temperature and metallicity in Fig. 1141 
are local quantities, and as such cannot be compared directly 
with those shown in Fig.[TS]in the following Section and with 
the average quantities inferred observationally for the ICM 
and IGM. This is important for the metallicity, as the metal 
distribution in our simulations is inhomogeneous on small 
scales. The scatter of the distributions, which is not shown, 
is relatively small at high fluxes and increases at lower fluxes 
to up to 0.5 dex. 

The trends observed in Fig. [14] are consistent with the 
results of Fig. 1131 The line flux is strongly correlated with 
density and metallicity, and the highest fluxes are always 
produced by dense and metal-rich gas near galaxies. These 
correlations are not surprising, because the emission scales 
with the element abundance and with the density squared. 
On the other hand, the median gas temperature that pro- 
duces a given flux is almost constant at the highest fluxes 
and quickly decreases at the lowest fluxes. This is a direct 
consequence of the fact that the line emissivity peaks in a 
narrow temperature range, as shown in Fig. [1] The specific 
value of the median temperature at the highest fluxes is al- 
most equivalent to the temperature where the line emissivity 
peaks, as also seen in Fig. 1131 This suggests that the detec- 
tion of a specific emission line yields the temperature of the 
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Figure 14. Median values of the gas particles's overdensity (left panel), temperature (middle panel) and metallicity (right panel) that 
produce a given X-ray particle flux at z = 0.25 for the lines CVI, OVII, OVIII and NeX. Note that fluxes are given in units of 
photon s~ 1 cm" 2 , but that previous figures used photon s _1 cm~ 2 sr" 1 . Because we plot the flux emitted per particle, it is directly 
dependent on the numerical resolution. The absolute normalization of the x therefore does not contain any useful information. For 
reference, p m ean(z = 0.25) corresponds to ripj ~ 4 X 10 _7 cm — 3 and our star formation threshold corresponds to log 10 p/p mcan Rj 5.4. 
Note that emission from higher density gas (i.e. the ISM) was ignored. The median density and metallicity increase with luminosity, 
while the median temperature flattens out at the peak temperature of the emissivity of each line. 



gas with a small degree of uncertainty. The value of the ob- 
served flux then provides information about the density and 
the metallicity of the gas. The degeneracy between the two 
might be hard to break without further data from multiple 
lines from different elements and ionisation states, but it ap- 
pears from Fig. [I4]that, at the bright end, the flux increases 
faster with the gas density than with the metallicity. 

For each line, the maximum flux is reached at the den- 
sity that corresponds to our star formation threshold. Be- 
cause we lack both the physics and the resolution to model 
the gas at higher densities (recall that we impose an effective 
equation of state onto the ISM), we have ignored emission 
from higher densities. The peak fluxes are therefore limited 
by this density cut. We stress that this means that the mod- 
els presented here underestimate the maximum fluxes. On 
the other hand, the higher fluxes arise in gas with densities 
typical of the ISM and can therefore hardly be counted as 
WHIM. 

Finally, we note that the trends with density and tem- 
perature remain the same if we assume a constant metal- 
licity rather than the metal distribution predicted by the 
simulation. 



6.3 Density and temperature cuts 

In this Section we create O VIII emission maps using only 
gas with temperatures and densities in a given interval. This 
will help us to understand the dependence of the emission 
on temperature and density and will give an indication of 
what type of gas produces what kind of emission. For com- 
parison, we show in Fig. [15] the density, temperature and 
metallicity of the gas in the same region of space shown in 
all the emission maps in this paper. 

Fig. [16] shows O VIII emission maps made using only 
gas with temperatures in a specific range. We consider the 
four temperature intervals T < 10 5 K, 10 5 K < T < 10 6 K, 
10 6 K < T < 10 7 K and T > 10 7 K. Similarly, in Fig. HTl we 
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Figure 18. As Fig. H but showing the PDF of the OVIII flux 
from gas in different temperature intervals. The highest O VIII 
flux is produced by gas with T > 10 6 K. O VIII emission is weak- 
est in the lowest temperature ranges. 



consider four different density intervals, with Log 10 p/p mean 
varying by one dex in each density bin. The flux PDFs for 
the entire simulation volume are shown in Fig. [18] for the 
temperature cuts and in Fig. \T§\ for the density cuts. Figs. 
U6H19l demonstrate that the strongest O VIII emission comes 
from gas in dense (Log 10 p/p mcan > 3) and hot (T > 10 6 K) 
regions. 
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Figure 15. As Fig. [3] but showing maps of the gas overdensity (left panel), temperature (middle panel) and metallicity (right panels). 
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Figure 16. As Fig. [3] but showing maps of O VIII emission from gas in four different temperature ranges. From left to right: T < 10 5 K, 
10 5 K < T < 10 6 K, 10 6 K < T < 10 7 K and T > 10 7 K. The highest O VIII fluxes are produced by gas with T > 10 6 K. O VIII emission 
is strongest in the highest temperature range. 



7 DEPENDENCE ON THE PHYSICAL 
MODELS 

A key feature of the OWLS project is the availability of a 
large number of runs with varying physical prescriptions. 
In this Section we investigate how a series of physical pre- 
scriptions affects the predictions for the WHIM emission. 
We will only show results for O VIII A18.97 A emission, but 
note that we find similar trends for other lines. We focus on a 
subset of the OWLS models with varying feedback schemes 
and cooling recipes, which are the physical implementations 
that produce the largest variations in the diffuse emission. 
These runs differ not only in the total amount of metals 
produced and in the spatial distr ibution of the met als, but 
also in the star formation history l|Schave et al.|[2O10l ) and in 



the density and temperature distribution of the diffuse gas. 
Physical prescriptions are changed one at a time, in order 
to isolate their impact on the results. The physical mod- 
els considered here are the following (we refer the reader to 
ISchave et all (|2010h and to the works cited below for more 
detailed descriptions) : 

(i) REF: Our reference model, as described in Section [2] 
All results presented in the preceding sections were ob- 
tained from this simul ation. This model uses the sta r for- 
mation prescription of ISchave fc Dalla Vecchial (l2008h. the 
supern ova-driven wind model of iDalla Vecchia fc Schavd 
(2008) with wind mass loading r\ = 2 and wind initial ve- 
lo city ti w = 600 kms" 1 , the metal-dependent cooling rates 
iersma et alJ l|2009a[ ) , the chemical evolution model de- 
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Figure 17. As Fig. [3] but showing maps of O VIII emission for gas in four different density ranges, as indicated at the top of each panel. 
The highest fluxes are associated with the densest gas. 
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Figure 19. As Fig. H but showing the PDF of the O VIII flux 
from gas in different density intervals. The highest fluxes are as- 
sociated with the densest gas. In particular, the flux predicted 
for gas outside the virial radii of haloes (i.e. p/pmcan < 200) is at 
least a few orders of magnitude weaker than the flux from denser 
regions. 



scribed in lWiersma et alJ l|2009bl ) and the WMAP3 cosmol- 
ogy. AGN feedback is not included. 

(ii) NOSN: As REF, but without supernova feedback. 
The metals produced by stars are only transported by gas 
mixing and no energy is transferred to the IGM and ISM 
when SNe explode. 

(iii) NOZCOOL: As REF, but with radiative cooling rates 
calculated assuming primordial abundances. The neglect of 
the contributions of metal lines implies slower gas cooling 
than in the REF model. 



(iv) AGN: As R EF, but with the additio n of AGN feed- 
back as described in iBooth fc Scharel \200$ ); 

(v) WML4 : As REF, but with a wind mass loading twice 
as high, that is r\ — 4. The total energy injected in winds is 
thus also twice that in REF; 

(vi) WDENS: As REF , but the wind prop erties scale with 
the local gas density as (|Schave et al.|[201QT ): 



; w = «w0 (p/Pcrit) 1 '' 6 
V = V0 (P/Pcrit)" 2/6 



(5) 
(6) 



where p cr i t is the density threshold for star formation, v v q = 
600 kms -1 and r]o = 2. This implies that at p cr i t the gas is 
kicked with the same velocity and mass-loading as in the 
REF model. This particular scaling of wind velocity with 
density was chosen so that the wind velocity scales with the 
local sound speed (recall that we impose the equation of 
state P oc p 4 / 3 onto the ISM) and the wind mass-loading is 
scaled such that the total wind energy is independent of the 
gas density. As a consequence, the wind energy is equal to 
that in the REF model. 

(vii) WVCIRC: As REF, but with "momentum- 
driven" galactic winds f ollowing the prescription of 
lOppenheimer fc Pavel ()2008l ). with the difference that winds 
are "local" to the star formation event and are fully 
hydrodynamically coupled a s in the REF model (see 
iDalla Vecchia fc Schave!l2008l for a discussion of the impor- 
tance of these effects). The wind initial velocity and the 
wind mass loading factor are defined as a function of the 
galaxy velocity dispersion a as v w = 5er and r\ = u w o/o" 
respectively, with w w o = 150 kms . The velocity disper- 
sion a — \/2«circ is estimated using an on-the-fly friends-of- 
friends halo finder. This model was motivated by the idea 
that galactic winds may be driven by radiation pressure on 
dust grains rather than supernovae. The energy injected into 
the wind becomes much higher than for REF for halo masses 
exceeding 10 11 - 10 12 M . 

(viii) MILL: As REF, bu t using the WMAP year 1 cos- 
mology dSpergel et alj|2003h as in the Millennium simulation 
l|Springel et al.ll2005h and a wind mass loading factor f] = 4. 
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Figure 20. As Fig. [3] but showing O VIII emission maps for different simulations, as indicated in the top left corner of each panel. 



Results from this model must be compared to those of the 
WML4 model to isolate the effect of cosmology. 

(ix) DBLIMF: As REF, but using a top-heavy IMF at 
high gas pressures. This model switches from a Chabrier 
IMF to a power-law IMF dN/dM oc M" 1 (as compared to 
•x M -2 - 3 for the high -mass tail of the Chabrier IMF) at 
the critical pressure P/k — 2.0 x 10 6 cm -3 K, which was 
chosen because for this value ~ 1CF 1 of the stellar mass 
forms at higher pressures. The extra SN energy is used to 
increase the initial wind velocity to 1618 kms -1 for star 
particles formed with a top-heavy IMF, while keeping the 
mass loading factor fixed at r\ = 2 (as for the Chabrier IMF, 
these parameter values correspond to 40% of the available 
SN energy). This model is called DBLIMFCONTSF VI 61 8 



in the OWLS project (which i ncludes more model s with top- 
heavy IMFs in starbursts, see lSchave et al.ll2O10T ). 

Fig. [20] shows maps of the O VIII emission for all runs 
in a region of 14 h^ 1 Mpc on a side or, equivalently, 1.12 
degrees on the sky. Fig. [2TJ shows the emission in the same 
simulations, but for a constant gas metallicity of Z = 0. IZq . 
Figs. [22] and [23] present the flux PDFs of the O VIII emis- 
sion using the predicted abundances and Z = O.IZq, re- 
spectively. 

The dominant impression we get from looking at these 
figures, is that most models predict very similar O VIII 
fluxes. This is somewhat surprising, as the physics that is 
included is quite different and the mo dels predict a rathe r 
wide range of star formation histories (Sc have et alj|201oh . 
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Figure 21. Same as Fig. 1201 but for constant gas metallicity Z = O.lZp 
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Closer inspection does, however, reveal several differences. 
Visually, we can easily see from the top panels of Fig. [20] 
that models NOSN, WVCIRC, and MILL differ most from 
our reference model. Comparing Figs. l20l and l21l we see that 
for NOSN the difference is driven by the metal distribution. 
Without galactic winds, the metals remain mostly confined 
to the inner parts of haloes and the flux from low-density gas 
is sharply reduced. The predicted metal distributions can, 
however, not account for the differences seen for WVCIRC, 
and MILL. The flux from high-density gas is considerably 
lower for WVCIRC, which indicates that the peak densities 
are lower and/or that the temperatures of high-density gas 
are lower than the value for which the emissivity of O VIII 
peaks (i.e. 10 6 3 K). The MILL run appears different mostly 



because structure formation has progressed further in this 
run, because it uses a higher value of erg. 

The flux PDFs shown in Figs. [22] and [23] confirm the 
conclusions drawn from the visual inspection of the emission 
maps. Looking first at the global distributions shown in the 
insets, we see that only NOSN differs substantially, the peak 
of the PDF having shifted to much lower fluxes. If a constant 
metallicity is assumed, then all models agree near the peak of 
the distribution and the low-flux tail of the PDF vanishes. 
This indicates that fluxes <C 10 -5 photon s _1 cm" 2 sr" 1 
imply gas metallicities <C O.IZq. For all models the peak 
fluxes are lower if a constant metallicity is assumed, from 
which we conclude that the high flux tail probes gas with 
metallicities in excess of 0.1/7© . 
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Figure 22. As Fig. [4] but showing the flux PDFs of O VIII emission for simulations with varying physical prescriptions. The simulations 
shown in the left panel are the same ones as in Fig. 1201 Additional runs are shown in the right panel and compared to the results of 
the REF run. The results are surprisingly robust to changes in the physical model. However, neglecting metal-line cooling results in a 
large overestimate of the maximum flux. Ignoring supernova-driven winds has little effect for the brightest pixels, but otherwise shifts 
the PDF to lower fluxes. 




Figure 23. As Fig. [4] but showing the flux PDFs of OVIII emission for simulations with varying physical prescriptions, under the 
assumption of a constant gas metallicity of Z = O.IZq. The simulations shown are the same ones as in Fig. 1201 The assumption 
of constant metallicity boosts the flux from the low-density, cool regions of the universe, where the largest contribution to the line 
emissivity comes from photo-ionisation by the UV background (compare the inset to that of Fig. 1221) . but it reduces the maximum fluxes. 
The differences between the models are qualitatively the same as when the predicted metal distributions are used (compare with Fig. 12211 , 
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The model for which the high-flux tail differs the most is 
NOZCOOL, which neglects metal-line cooling. Its flux PDF 
starts to differ from that of REF for fluxes > 1 photon s _1 
cm" 2 sr" 1 , reaching a maximum flux that is about an or- 
der of magnitude higher. The fact that metal-line cooling 
is important is not surprising, as the metal-line emission 
that we are studying dominates the radiative cooling rates 
for temperatures ty pical of the WHIM a nd for metallici- 
ties > 10" 1 Z Q (e.g. IWiersma et al.ll2009ah . Ignoring metal- 
line cooling during the simulations leaves the gas in the 
temperature-density regime where it emits efficiently, artifi- 
cially boosting the predicted fluxes. The calculation of the 
emission from metals in this model is of course inconsistent 
with the assumption of primordial abundances used to cal- 
culate the gas cooling rates. However, this in consistency is 
present in other work (e.g. iFang et ail 120051 ) . which often 
used cooling rates based on primordial abundances during 
the simulations and calculates the metal-line emission of the 
gas in post-processing. Clearly, this approach will strongly 
overestimate the maximum fluxes. 

While the emission map of the MILL model appeared 
to be substantially different (Fig. [20}, the change in cosmol- 
ogy does not change the statistics of the O VIII emission 
substantially, as can be seen by comparing the flux PDFs of 
the WML4 and MILL models. 

There are enormous differences in the feedback pro- 
cesses that are included in the simulations shown here, rang- 
ing from no galactic winds (NOSN), to efficient winds driven 
by supernovae ( WDENS), to extremely strong winds in star- 
bursts resulting from a top-heavy IMF (DBLIMF), to vio- 
lent AGN feedback (AGN). These differences s trongly affect 
the star formation histories |Schave et al.|[201ol ). yet the flux 
PDFs of all these models are strikingly similar. This sug- 
gests that the metal distributions are also similar in the gas 
that emits O VIII efficiently. Indeed, we will show elsewhere 
(Wiersma et al., in preparation; McCarthy et al., in prepara- 
tion) that the metallicity of the WHIM (much more so than 
that of other components) is remarkably robust to changes 
in the feedback prescription. Indeed, assuming a constant 
metallicity does not change the relative difference between 
the models much. This suggests that the differences in the 
high flux tails are driven by differences in the distribution of 
the baryons in the temperature (T ~ 10 6 ' 5 K) and density 
(p 2> 10 3 p mC an) regime that dominates the emission from the 
brightest regions (see Fig. 1141) . In particular, it appears that 
the "momentum-driven" winds ( WVCIRC), which predict 
substantially lower peak fluxes (even for a constant metal- 
licity), substantially reduce the peak gas densities of the 
WHIM. This suggests that these winds are able to remove 
more gas from the centers of haloes, perhaps because the 
winds in this model use much more energy than is available 
from supernovae. 

Thus, we conclude that bright O VIII emission, and 
other soft X-ray lines too, may tell us more about the gas 
densities in the centers of haloes than about the metallicities 
and the temperatures (which must be close to the peak of 
the emissivity curve for the emission to be bright). This is 
because the emissivity scales as Zp 2 and is thus most sensi- 
tive to the gas density. 



8 CAN THE WHIM BE DETECTED IN 
EMISSION? 

In this Section we discuss the possibility that current and 
future X-ray telescopes may detect the WHIM emission at 
the levels predicted by our simulations. 

Previous studies predict that X-ray emission from the 
WHIM is far below the detection threshold of the Chan- 
dra and XMM-Newton telescopes (lYoshikawa et ail 120031 ; 
lYoshikawa et all |2004| ; IFang et ail 120051 ). Predictions that 
WHIM emission lines in filaments co uld be detected b y 
XMM-Newton have not been confirmed dPierre et ail 2000). 
However, future mission s such as EDGE ( Piro et ail 12003 ). 



IXO and Generation-X ^Windhors t et ail 



20061 ). have been 



shown to have a better chance of detect i ng diffuse emiss ion 
l|Yoshikawa et aill2004 IFang et al.ll2005l ; lBregmanll2007l ). 

Ideally, the emission from the WHIM requires an in- 
strument which couples a large field of view (FoV, hereafter) 



with high angular and s pectral resolution ( Yoshikawa et ail 
|2004 IFang et ail 120051 ; IPaerels et ail 12008 ). While instru- 



ments with high angular and spectral resolutions alone are 
fine for searches of absorption features, a large FoV is desir- 
able for mapping diffuse emission. High angular resolution is 
also very important, as we found that the emission is domi- 
nated by bright, compact regions. Without sufficient angu- 
lar resolution, it will not only be more difficult to detect the 
emission, but one would also risk misinterpreting its origin. 
For example, emission that looks smooth and filamentary at 
poor angular resolution may come from unresolved, dense 
knots rather than from diffuse, low density gas. 

A few dedicated small to medium-size missions have 
been proposed to investigate the WHIM emiss ion. The Miss- 
ing Baryon Explorer (MBE, IFang et aill2005l ) was the first 
of such medium-size missions, later fol lowed by the Diffus e 
Intergalactic Oxygen Surveyor (DIPS, [Ohashi et ail [2006). 
the Extreme phySics in the TR ansient and Evolving cos MOs 
""' 20061) . EDGE (|Piro et aill2009l) and 
2009). The last two missions are 



( ESTR EMO,\Piro et al 



Xenia (|Hartmann et al 



very similar. The common goal of all these missions is to 
build an instrument with a large FoV and high spatial and 
spectral resolution. 

The proposed technical specifications of the instruments 
we consider are listed in Table [3] The FoV of the X-ray 
Microcalorimeter SpectrometeiH (XMS) on board IXO is 
5.4' x 5.4', which corresponds to 1.1 ft -1 comoving Mpc 
at z = 0.25. Als o on IXO, the Wide Field Imager (WFI, 
iTreis et alj 12009 ) has a FoV of 18' x 18', corresponding to 
3.7 h^ 1 comoving Mpc at the same redshift. Both instru- 
ments have sufficient coverage to image gas in filaments at 
z ~ 0.25, with an angular resolution of 5" each. For com- 
parison, the Wide Field Imager (WFI) and the Wide Field 
Spectrometer (WFS) on Xenia were proposed to have FoVs 
of 1 and 0.7 degrees and angular resolutions of 15" and 3.7', 
respectively. While XMS on IXO and WFS on Xenia have 
impressive energy resolutions of about 3 eV, it is much worse 
for the WFIs (150 eV for IXO and 70 eV for Xenia). 

The unresolved cosmic X-ray background has been mea- 
sured to be ~ 13 ± 1 photon s -1 cm" 2 sr - 1 keV -1 in 
the range 0.5 - 1 keV l|Hickox fc Markevitchl 120061 ) and 
is thought to be dominated by diffuse Galactic and local 



http: / /ixo. gsfc.nasa.gov/technology/xms. html 
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Table 3. Summary of the technical characteristics of the instruments on IXO and Xenia. 



Instrument 


Effective area 


FoV 


Angular resolution 


Energy range 


Energy resolution 




(m 2 ) 


(arcmin 2 ) 


(arcsec) 


(keV) 


(eV) 


IXO 












XMS 


3 © 1.25 keV 


5.4 X 5.4 


5 


0.3 - 7.0 


2.5 @ central 2' x 2' 


WFI 


3 @ 1.25 keV 


18 X 18 


5 


0.1 - 15 


150 


Xenia 












WFS 


0.12 @ 0.6 keV 


42 X 42 


222 


0.1 - 2.2 


3 @ 0.5 keV 


WFI 


0.058 ® 1 keV 


60 x 60 


15 


0.3 - 6 


70 & 1 keV 



thermal-like emission. For the energy resolution of the XMS 
(AE = 2.5 eV) this corresponds to about 3 x 10 -2 photon 
s _1 cm~ 2 sr —1 . 

The instrumental background for XMS is about 2 x 
10 -2 photon s _1 keV -1 per cm -2 of the detector area (Kaas- 
tra 2010, private communication). For a pixel size of 300 fj,m 
and 3 arcsec this corresponds to 8.5 x 10 11 photon s _1 sr -1 
keV -1 . For an effective area of 2 x 10 4 cm 2 and an energy 
resolution of 2.5 eV this implies a limiting surface bright- 
ness of 1.1 x 10~ 2 photon s" 1 cm 2 sr 1 , slightly below the 
extragalactic background. Thus, we expect observations to 
become background limited below 10" 1 photon s" 1 cm~ 2 
sr -1 . For an angular resolution of 15 arcsec and an effective 
area of 2 x 10 4 cm 2 , this surface brightness would result in 
about 10 detected photons per resolution element in a 1 Ms 
exposure. Hence, with an instrument like XMS it would be 
possible to detect surface brightnesses > 1 photon s _1 cm -2 
sr -1 (at 15 arcsec resolution) with high confidence in reason- 
able exposure times. Fig.[l9]shows that this would enable the 
detection of O VIII from gas with p > 10 2 p mC an at z = 0.25. 

For the proposed WFS for Xenia the instrumental back- 
ground is expected to be even lower and the energy resolu- 
tion similar to that for XMS on IXO. The angular resolution 
is, however, much poorer (about 4') and the effective area is 
an order of magnitude smaller (10 3 cm 2 ). A 1 Ms exposure 
would in that case yield ~ 10 2 photons per resolution ele- 
ment for a surface brightness of 10 _1 photon s _1 cm -2 sr -1 , 
more than sufficient for a robust detection. Note that XMS 
on IXO would detect an order of magnitude more photons 
when averaged over the same angular size (because its effec- 
tive area is an order of magnitude larger). We predict from 
the equivalent version of Fig. [19] for a pixel size of 4 arcmin 
(not shown) that this surface brightness limit is sufficient to 
detect some O VIII emission from gas with p > 10 3 p mC an at 
2 = 0.25. Note that while this emission would appear to be 
diffuse due to the low angular resolution, it would in fact 
come from unresolved, dense clouds. 

Compared with the WFS, the WFI for Xenia has much 
better angular resolution (15"), a smaller effective area 
(about 0.6 x 10 3 cm 2 ) and a much worse energy resolution 
(AE — 70 eV). Because the spectral resolution is more than 
an order of magnitude lower than for the WFS, the back- 
grounds will be higher by the same factor. For reasonable 
exposure times (~ 1 Ms) we would still not be background 
limited. The poor spectral resolution would, however, make 
it much more challenging to interpret the observations, as 
the lines correspond to the differential Hubble flow across a 
comoving distance ~ 250 hT 1 Mpc (for z = 0.25). On the 



other hand, the interpretation would be greatly facilitated 
by the much improved angular resolution, although this will 
not necessarily help with the identification of the lines. 

The WFI for IXO has an energy resolution that is lower 
still, AE = 150 eV, which may lead to severe projection ef- 
fects and which would make it challenging to identify the 
lines. While WFI does not have a better angular resolution 
than XMS, it does have a much larger field of view (18 ar- 
cmin vs. 2 arcmin for the high energy resolution mode of 
XMS). 

According to our predictions, the WFS on Xenia should 
be able to detect strong emission from the WHIM only for 
gas with overdensities greater than a thousand, such as in 
the centres of groups. IXO would be extremely exciting as 
it would be sufficiently sensitive to detect emission from 
O VIII, and probably a few other lines, from both the cores 
of small groups and from the outskirts of larger groups and 
clusters down to overdensities > 10 2 , and with an angular 
resolution that is sufficient to resolve much of the intrinsic 
structure. 

Instruments with higher sensitivity and lower flux limits 
will be required in the future to image the emission from the 
diffuse WHIM in low-density regions and in filaments where 
p <C 10 2 p mC an- Generation-X (lWindhorst et al.ll200a l might 
be one of those missions. Since most of the baryons contained 
in the WHIM reside at such low overdensities (Fig. I12|l , it 
seems unlikely that the bulk of the baryons will be detected 
in soft X-ray emission any time soon. Absorption experi- 
ments clearly have a better chance to accomplish this goal. 

However, we have seen that detecting soft X-ray emis- 
sion from the gas in and around (groups of) galaxies is within 
reach of proposed missions. This is extremely exciting, as 
such observations would provide vital clues on gas infall, 
outflows and metal enrichment, which are the key, poorly 
understood processes in models of the formation of galaxies. 



9 CONCLUSIONS 

Gravitational accretion shocks and, to a lesser extent, galac- 
tic winds driven by supernovae and/or AGN, can shock-heat 
the diffuse intergalactic medium to temperatures > 10 J K. 
At low redshifts a large fraction of the baryons are thought 
to reside in this warm-hot intergalactic medium (WHIM; 
10 s K < T < 10 7 K) and soft X-ray emission lines from met- 
als offer a way to detect them. 

Here, we used a subset of cosmological simulations from 
the Overwhelming ly Large Simulations (OWLS) project 
ijSchave et al.l[2oTol ) to study the nature and detectability of 
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this emission. These simulations are among the largest dissi- 
pative runs ever performed, and use new prescriptions for a 
range of baryonic processes, such as star formation, galactic 
winds driven by massive stars and supernovae, and feedback 
from AGN. Of particular relevance for us, is that they are 
fully chemodynamical, tracking the delayed release of 11 dif- 
ferent elements by massive stars, supernovae of types II and 
la, and asymptotic giant branch stars, and that radiative 
cooling is implemented element-by-element in the presence 
of an evolving UV/X-ray background. 

Our main conclusions can be summarised as follows: 

(i) The strongest emission in the soft X-ray band comes 
from the O VIII line, whose emissivity peaks at a tempera- 
ture of a few times 10 6 K. A large number of other lines are, 
however, only slightly weaker. 

(ii) Emission lines from metals provide strong constraints 
on the temperature of the emitting gas. Most of the emit- 
ted radiation and essentially all of the potentially detectable 
flux comes from gas with temperatures close to the peak of 
the emissivity curve for a collisionally ionised gas. Because 
of the quadratic dependence of the emissivity on density, 
the emission traces neither the bulk of the baryonic nor the 
bulk of the metal mass. Instead, the emission is dominated 
by metal-rich gas at very high overdensities in and around 
galaxies. Since the peak temperature of the emissivity in- 
creases with atomic number and ionisation state, lines from 
heavier or more highly ionised atoms trace hotter gas. As 
such, C VI and O VII trace gas with T ~ 10 6 K, O VIII, 
NelX and FeXVII gas with 10 6 K < T < 10 7 K, and NeX, 
MgXII, Si XIII and S XV gas with T ~ 10 7 K. 

(iii) By comparing results from different OWLS models, 
we have investigated the dependence of the results on a 
number of physical mechanisms. It is essential to accurately 
model the cooling rates. In particular, studies that neglected 
metal-line cooling are not self-consistent and will have over- 
estimated the maximum fluxes by about an order of mag- 
nitude. Without galactic winds, the PDF of the flux shifts 
to much lower values. However, the high-flux tail, which will 
become detectable first, is surprisingly robust to strong vari- 
ations in the strength of the winds, remaining similar even 
when galactic winds are turned off or when violent AGN 
feedback is included. While assuming a constant metallic- 
ity of O.IZq reduces the width of the flux PDF (and thus 
suppresses the maximum fluxes), the relative differences be- 
tween the models remain similar. This suggests that the dif- 
ferent feedback processes mainly affect the metal-line emis- 
sion by changing the gas densities in and around galaxies, 
as the highest fluxes always come from regions with tem- 
peratures close to the peak of the emissivity curve and with 
(local) metallicities > Zq. 

(iv) Soft X-ray emission lines in the low-redshift Universe 
are likely to be detectable only for relatively high-density re- 
gions, such as groups and the outskirts of clusters. Because 
dense, compact structures can easily dominate the emission 
even if they contain only a very small amount of mass, it 
is important for future observatories to have high angular 
resolution. Without sufficient angular resolution, the emis- 
sion is harder to detect and one risks misinterpreting the 
origin of the detected emission. For example, emission that 
looks diffuse may in reality come from a number of compact 



gas clouds in and around galaxies that trace the large-scale 
structure. 

(v) While soft X-ray line emission is too faint for it to be 
a promising route to close the baryon budget, it does offer 
the exciting possibility to image the gas accreting onto and 
flowing out of galaxies. Currently proposed missions could 
provide vital clues on some of the most poorly understood, 
but crucial ingredients of models of the formation and evo- 
lution of galaxies. 
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APPENDIX A: CONVERGENCE TESTS 

In this Appendix we present a number of tests we have per- 
formed to verify how numerical parameters affect our re- 
sults. In particular, we consider the effect of resolution in 
Appendix IA1I Results for different box sizes (at constant 
resolution) are shown in Appendix I A2I In Appendix [A3] we 
investigate the effect of the slice thickness on the shape of 
the flux PDF. All tests use the REF physical model and 
results are shown only for the O VIII 18.97 A line. 



Al Resolution 

To test how the numerical resolution affects our results, 
we have created O VIII emission maps for three differ- 
ent realisations of the same simulations with particle num- 
bers N = 2 x 128 3 , 2 x 256 3 and the default number of 
2 x 512 3 . The corresponding dark matter particle masses 
are 2.6 x 10 10 , 3.2 x 10 9 and 4.1 x 10 s h' 1 M , respectively. 
The corresponding maximum physical gravitational soften- 
ing lengths are 8, 4, and 2h _1 kpc. All runs use a box size 
of 100 h" 1 comoving Mpc. 
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Figure Al. As Fig. [4] but showing the O VIII flux PDF as a 
function of the resolution of the simulation. The particle mass 
(gravitational softening) increases by a factor of eight (two) when 
the number of particles decreases by a factor of 2 3 . All runs as- 
sume the same physical model, angular resolution and box size. 
The predictions of the N = 256 3 and 512 3 runs are very close for 
fluxes greater than 10 -5 photon s _1 cm -2 sr -1 . 

Fig. IA1I shows the flux PDF for maps with 15" angu- 
lar resolution. The lowest resolution run with N = 2 x 128 3 
deviates most from the other two, in particular for fluxes 
lower than about 1 photon s _1 cm -2 sr -1 . This is at least 
in part caused by the fa ct that this run se verely underes- 
timates the stellar mass <|Schave et al.l l2010) and hence the 
gas metallicity. The higher resolution runs with 2 x 512 3 and 
2 x 256 3 particles show good agreement, although there are 
small differences at the highest fluxes. 



A2 Box size 

We have tested the effect of the size of the simulation box on 
the flux PDF of O VIII emission, while keeping the resolu- 
tion constant. Fig. IA2I compares simulations implementing 
the REF model with box sizes of 100, 50 and 25 h' 1 Mpc, 
respectively, and with constant dark matter particle mass 
772dm =4.1 X 10 8 ft" 1 Mq. The flux PDFs converge for in- 
termediate flux values, but diverge at the lowest fluxes. The 
lowest and very highest fluxes are only found in the largest 
box, which contains the greatest number of particles and 
provides the most pixels. It thus samples the tails of the 
distribution best. 



A3 Slice thickness 

In this Appendix we describe how varying the thickness of 
the slices in which we divide the simulation box affects the 
flux PDF. In Fig. IA3I we show this by comparing the flux 
PDFs, normalised by the slice thickness, for the four differ- 



Figure A2. As Fig. g] but showing the OVIII flux PDF for 
simulations with box sizes of 100, 50 and 25 h~ 1 Mpc. All runs 
assume the same physical model, particle mass and angular res- 
olution. Except for very low fluxes (< 10 — 5 photon s — 1 cm -2 
sr -1 ), the results are insensitive to the size of the simulation box. 
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Figure A3. The O VIII flux PDFs, normalised by the slice thick- 
ness, for maps with varying slice thickness: 40, 20, 10 and 5 h~ 1 
Mpc. The number of slices used to calculate each PDF is 2, 5, 
10 and 20, respectively. The test is for the REF model with box 
size 100 /i -1 Mpc at z = 0.25. All maps assume the same angular 
resolution of 15" . In the high flux regime the PDF is proportional 
to the slice thickness. 
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ent, evenly-spaced slice thicknesses of 40, 20, 10 and 5 h~ x 
Mpc. 

In the potentially detectable regime, the results appear 
to be independent of the slice thickness when the flux PDF 
is normalised by the slice thickness. In other words, the PDF 
in the high flux regime is proportional to the slice thickness. 
This is expected if the slice is sufficiently thick to fully con- 
tain the objects giving rise to the emission, but sufficiently 
thin for the superposition of multiple structures to be unim- 
portant. 

The shape of the PDF at low fluxes, instead, depends on 
the level of the flux in low-density regions and on the num- 
ber of low-density structures such as filaments that could be 
found in each slice, which is dependent on the slice thick- 
ness. As a consequence, the actual width of the PDF at 
intermediate fluxes scales with the inverse of the slice thick- 
ness and the number of pixels at the lowest fluxes increases 
for decreasing slice thickness. 

When creating the slices, we neglect the peculiar ve- 
locities of particle, selecting them only according to their 
position within the simulation box. This is justified for a 
large slice thickness, when the peculiar velocities are small 
compared with the difference in the Hubble flow across the 
slice, because the effect of particle displacement is small. 
However, it may have a small effect on the predicted fluxes 
for the 5 /i -1 Mpc slices, which corresponds to 443 kms -1 
at z = 0.25. 

Throughout this work, we have assumed a slice thick- 
ness of 20 /i -1 Mpc, but Fig. IA3I demonstrates that the re- 
sults can easily be scaled to other slice thicknesses. 

This paper has been typeset from a TjnjX/ DTj^X file prepared 
by the author. 
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